clc;clear;close all;
B = linspace(50,120,8);
N = linspace(3,10,8);
F = linspace(10,80,8);
X = [B;N;F];
strX = ['B','N','F'];
RMSE = zeros(4,8);
TIME = zeros(6,8);
for i = 1:size(X,1)-2
    for j = 1:size(X,2)
        path = ['Data\',strX(i),'\',strX(i),'=',int2str(X(i,j)),'\LS\'];
        rmse_LS = load([path, 'rmse.mat']);RMSE(1,j)  = rmse_LS.rmse;
        time_LS = load([path, 'time.mat']);TIME(1,j) = time_LS.time;
        
        path = ['Data\',strX(i),'\',strX(i),'=',int2str(X(i,j)),'\WLS\'];
        rmse_WLS = load([path, 'rmse.mat']);RMSE(2,j)  = rmse_WLS.rmse;
        time_WLS = load([path, 'time.mat']);TIME(2,j) = time_WLS.time;
        noiseTime_WLS = load([path, 'noiseTime.mat']);TIME(3,j) = noiseTime_WLS.Time;
        
        path = ['Data\',strX(i),'\',strX(i),'=',int2str(X(i,j)),'\WLS(Zhang)\'];
        rmse_WLSZhang = load([path, 'rmse.mat']);RMSE(3,j)  = rmse_WLSZhang.rmse;
        time_WLSZhang = load([path, 'time.mat']);TIME(4,j) = time_WLSZhang.time;
        noiseTime_WLSZhang = load([path, 'noiseTime.mat']);TIME(5,j) = noiseTime_WLSZhang.Time;

        path = ['Data\',strX(i),'\',strX(i),'=',int2str(X(i,j)),'\WPCA\'];
        rmse_WPCA = load([path, 'rmse.mat']);RMSE(4,j)  = rmse_WPCA.rmse;
        time_WPCA = load([path, 'time.mat']);TIME(6,j) = time_WPCA.time;

    end
end

plot(B,RMSE(1,:), B,RMSE(2,:),B,RMSE(3,:),B,RMSE(4,:));

% plot(B,TIME(1,:), B,TIME(2,:),B,TIME(3,:),B,TIME(2,:)+TIME(3,:));